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Abstract: We show that mass varying neutrinos (MaVaNs) can behave as a negative pres- 
sure fluid which could be the origin of the cosmic acceleration. We derive a model independent 
relation between the neutrino mass and the equation of state parameter of the neutrino dark 
energy, which is applicable for general theories of mass varying particles. The neutrino mass 
depends on the local neutrino density and the observed neutrino mass can exceed the cosmo- 
logical bound on a constant neutrino mass. We discuss microscopic realizations of the MaVaN 
acceleration scenario, which involve a sterile neutrino. We consider naturalness constraints 
for mass varying particles, and find that both eV cutoffs and eV mass particles are needed to 
avoid fine-tuning. These considerations give a (current) mass of order an eV for the sterile 
neutrino in microscopic realizations, which could be detectable at MiniBooNE. Because the 
sterile neutrino was much heavier at earlier times, constraints from big bang nucleosynthesis 
on additional states are not problematic. We consider regions of high neutrino density and 
find that the most likely place today to find neutrino masses which are significantly different 
from the neutrino masses in our solar system is in a supernova. The possibility of different 
neutrino mass in different regions of the galaxy and the local group could be significant for 
Z-burst models of ultra-high energy cosmic rays. We also consider the cosmology of and the 
constraints on the "acceleron" , the scalar field which is responsible for the varying neutrino 
mass, and briefly discuss neutrino density dependent variations in other constants, such as 
the flne structure constant. 
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1. Introduction 

In recent years, we have learned much about the composition of the universe. Motivated ini- 
tially by the evidence from supernovae that the universe is accelerating [1,2], and with recent 
data on the cosmic microwave background radiation (CMBR) [3-5], the ACDM cosmology 
has become firmly established as the benchmark against which all other cosmologies must 
be compared. The need for a mysterious new "dark energy" component to the universe has 
pointed to the existence of new physics beyond the standard model (BSM). 

At almost the same time, the field of neutrino masses has undergone a revolution. There is 
the compelling evidence for atmospheric neutrino oscillations from the SuperKamiokande ex- 
periment [6], consistent with earlier evidence from the 1MB experiment [7,8] and Kamiokande 
[9,10], and recently confirmed by K2K [11-13]. The SNO [14] and KamLAND [15] exper- 
iments have confirmed the neutrino oscillation interpretation of the solar neutrino deficit 
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observed at other experiments [16-20]. The LSND result [21] remains a puzzling anomaly, 
whose interpretation via neutrino oscillations will be tested at the upcoming MiniBooNE ex- 
periment [22]. To add neutrino masses into the standard model, one must invoke either the 
existence of non-renormalizable operators, or light, sterile (singlets under the standard model 
gauge group) fermions, cither of which point to the existence of BSM physics. 

While both of these developments arc tremendously exciting, giving the second and third 
direct pieces of evidence for new physics (the first being the large body of evidence for non- 
baryonic dark matter), they are both unsatisfying. Models of dark energy are difficult to 
impossible to test. Some give precisely a small cosmological constant as a result, others 
predict a non-standard equation of state, but little else to distinguish them. A notable 
exception has measurable long-distance modifications of gravity as a consequence [23], but a 
theoretical debate still rages around this model [24,25]. Effective field theory considerations 
motivate modifications of gravity at the sub-millimeter scale [26,27], but it is not obvious 
that we will be able to discern the nature of the dark energy, even in the presence of such 
measurements. 

The situation for neutrino masses is not much better. Although the existence of neutrino 
mass is very exciting, it is not a priori clear how much we can learn about fundamental physics 
from the neutrino masses. Majorana neutrino masses can arise from a dimension five operator, 
whose origin can be from physics at an inaccessibly high scale, as high as lO^^GeV. Although 
there are exciting ideas that relate neutrino masses to supersymmetry breaking [28-32] , or the 
size of large extra dimensions [33,34], the most popular seesaw models [35,36] seem impossible 
to test directly. 

1.1 Coincidences and Damn Coincidences 

The problem of the dark energy is especially puzzling because of what has become known as 
the "cosmic coincidence problem." Namely, at the present time, pcdm/pa ~ 1/3, although 
this ratio changes as a function of cosmic scale factor as Hence, although these energy 

densities vary in dramatically different fashion over the history of the universe, we find that 
they are the same order at the present day. 

This coincidence has been the focus of a great deal of work: namely how to arrange for 
the dark matter and dark energy densities to be comparable. Solutions to this are difficult in 
that we actually know a great deal about the cosmological behavior of dark matter. Results 
from CMB let us know that at z ~ 1100 the universe was dominated by dark matter, and 
given the observed density of dark matter, this is consistent with a fluid rcdshifting with 
equation of state w = 0. While we know little about the dark energy at early times, we 
know that it is presently redshifting with equation of state w ^ —1 [37]. Thus it either has 
had u; — 1 for a long (cosmological) time, in which case there is a coincidence problem, or 
it tracked dark matter and only recently switched to redshift slowly, in which case one has 
a "why-now?" problem. In a sense, the cosmic coincidence between dark matter and dark 
energy might be more properly phrased as to why the dark energy now should be precisely 
three hundred million times smaller than the dark matter energy density at recombination. 
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Figure 1: Cosmological energy densities as a function of redshift in the ACDM model. The different 
lines are: long dashed - CDM; long-short dashed - baryons; short-dashed - neutrinos; solid sloped - 
radiation; solid horizontal - cosmological constant. 

On top of this, we are in fact confronted with a number of coincidences. There are many 
different components of the energy density of the universe: dark matter, baryons, radiation, 
neutrinos as well as dark energy. The first four all redshift very differently from with the 
dark energy, and yet we have the remarkable coincidence that they all have been equal to 
the dark energy within a redshift of a few, shown in figure |l]. Is it possible that another of 
these other coincidences is not a coincidence? While dark energy cannot reasonably track 
dark matter since recombination, is it possible that dark energy has tracked one of these 
other components? We have already discussed the difficulties in explaining the dark matter 
coincidence. The baryon number of the universe, constrained considerably by WMAP, is 
already probed at an even earlier era via the success of big bang nucleosynthesis (BBN), 
meaning that dark energy should not have been tracking baryons, unless it has the "why- 
now?" late time transition of its equation of state. Photon radiation is just as problematic. 

Neutrinos, however, are an entirely different story. Currently, the neutrino energy density 
is quite uncertain. Large scale surveys and WMAP together constrain ^l^y ^ .02. Terrestrial 
measurements of neutrino mass indicate ^1^ > 7 x 10~^. BBN implies that roughly three 
species of neutrino were relativistic at BBN, which leads to the mild conclusion that neutrino 
masses are below O(MeV). In reality, very little is known about the cosmological behavior of 
neutrinos and the neutrino energy density. 

In this paper, we will propose that the coincidence ~ p\ (within a factor of 10^) is 
not a coincidence at all, but in fact a relationship which holds over a large portion of the 
history of the universe. This is possible if the mass of the neutrino is not a fixed parameter, 
but a dynamical quantity, with an associated potential, in analog with variable mass dark 
matter [38]. The potential will have a zero energy minimum at a value of mu which is 
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larger than its present value. However, the presence of a (cosmologically) uniform neutrino 
background will lead to an effective potential which prevents rrii, from becoming too large, 
leaving a homogeneous negative pressure fluid in the universe (i.e., dark energy). 

Variable mass particles have been previously considered in many contexts, with attention 
recently revived in light of the current cosmological data. Early consideration of such ideas for 
dark matter with varying mass appeared in refs. [39,40] in the context of scalar-tensor gravity. 
Variable mass particles were proposed to solve the age of the universe problem in ref. [38], and 
more recently as a solution to the problem of the coincidence of dark matter and dark energy 
densities [41]. Phenomenological problems with this have been studied in [42-44]. In the 
context of neutrinos, ref. [45] considered a Yukawa- type coupling to an extremely light scalar 
field, leading to a different evolution of structure, and noted that the neutrino mass would 
be density dependent, although heavier rather than lighter at earlier times. Ref. [46] argued 
that dark energy in a desirable range could be linked to a scalar field coupled to neutrinos 
with a desirable mass. Refs. [47, 48] considered a neutrino Yukawa interaction as a source 
for "neutrino clouds" which might seed star formation, and change the mass of neutrinos 
in their vicinity. Ref. [49] considered coupling a sterile neutrino to a slowly rolling, Hubble 
mass scalar field which was responsible for dark energy, but did not consider the impact of 
the cosmic neutrino background on the potential. Ref. [50] considered that varying neutrino 
mass might have an impact on electroweak baryogenesis. 

In this paper we consider the possibility that the neutrino mass arises from an interaction 
with a scalar field, the "acceleron", whose effective potential changes as a function of the 
neutrino density, and consider the conditions under which such a new "dark force", felt by 
neutrinos, can explain the observed acceleration of the universe. In section ^ we will introduce 
the formalism and general features for mass varying neutrinos (MaVaNs), focusing on model 
independent relationships. We find a relationship between the neutrino mass and density, 
the neutrino dark energy density, and the equation of state of neutrino dark energy which 
will hold for any theory of neutrinos of varying mass. We discuss the effects of varying 
neutrino mass, specifying where we expect significant variations both in time and space, and 
reconsidering the traditional neutrino mass constraints as they apply to variable masses. We 
find that some neutrino mass bounds, such as those from large scale structure formation, do 
not apply. In section ^ we consider the details of particular microscopic realizations. Radiative 
stability of the scalar potential uniquely chooses theories with light particles and a cutoff for 
radiative effects at low energies (i.e., sub-eV). We argue that variable mass WIMPs are 
too heavy to be significant contributors to the dark energy, phenomenological considerations 
aside. A simple realization of MaVaNs arises from integrating out a standard model singlet 
fermion whose mass depends on the expectation value of the acceleron field. We note that 
traditional cosmological constraints from BBN on sterile neutrinos do not apply here, and 
that naturalness suggests that MiniBooNE will be able to detect such a sterile neutrino. In 
section ^ we consider phenomenological implications of MaVaNs, including effects at terrestrial 
neutrino experiments, as well as in solar oscillations and supernovae. 
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2. MaVaNs 



MaVaNs have a number of significant phenomenological consequences, and a great deal of this 
scenario can be understood in a completely model independent way. To begin, wc consider 
the neutrino mass nii, to be a dynamical field. Note that although for some purposes one 
should study m,^ as a function of a canonically normalized scalar field A, as long as dmn/dA 
does not vanish, for the purpose of analyzing the minimal energy density it suffices to write 
the scalar potential as a function of m^. For simplicity, we only consider a single nonvanishing 
neutrino mass, although generalization to more than one mass is straightforward. 

In general, the contribution of a neutrino background to the energy density is given by 



/rjit. , 
-^^W^Jik) , (2.1) 



where f{k) is the sum of the neutrino and anti neutrino occupation numbers for momentum 
k. 

The dependence of this energy density on a change in the neutrino mass is given by 

"'"'^s.^f^^^m- (2.2) 



Sm„ J {27r)3 ^e+mf 

The quantity , which we refer to as the "scalar neutrino density" , turns out to be important 
for MaVaN cosmology and phenomenology as the source term for a scalar field is proportional 
to Sjy. In the nonrclativistic limit s,y = n,y. the total density of neutrinos and anti neutrinos. 

For the cosmological Big Bang remnant neutrinos we have a contribution to the effective 
potential 

sv= ! ^TFTT^g 

J (27r)3i + exp(/3ov^^4TfeV)' 

where [3q = I/ZcTq is the temperature at freeze out, rriQ is the neutrino mass at freeze out, and 
z is the redshift at freeze out. In the nonrclativistic limit, the leading mu dependent term in 
the energy density from the neutrino background is simply n^rriiy. Therefore, in this limit the 
dependence of the energy density on rrii, is given by the effective potential 

V{m„) = m„n„ + Vb(m,,) , (2.4) 

where Vb(miy) is the scalar potential. The presence of thermal background neutrinos acts as 

a source which will drive rriy to small values. On the other hand, we assume Vq is minimized 
for a large value of m^^. Thus, these two terms compete, with a minimum at an intermediate 
value of 777-1/, with non-zero value for Vq. As the universe expands, the neutrino background 
dilutes, decreasing the source term, and 777^ is driven to larger values. 
At the minimum of the effective potential, 

V'{my) = n^ + V^{my) = . (2.5) 
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The equation of state of the total energy density from this sector is 



V + 
V ■ 



(2.7) 
(2.8) 



These equations already contain a great deal of information. Since the observed equation 
of state is close to w = —1, we know that the energy density in neutrinos alone should be 
small compared to the energy density in the total dark energy sector. This corresponds to a 
fairly flat potential Vo(m,y), either because the potential for A is rather flat, or because the 
dependence of rrii, on A is steep. Note however, that the constraints of the potential are not 
nearly as severe as, e.g., those in models of quintessence, where the ratios of derivatives of 
the potential to the potential itself are required to be smaller than the Hubble scale. Good 
candidate potentials Vo{mi,) grow as small fractional powers or logarithms of the neutrino 
mass. 

To the extent that dw/dz is small, or, equivalently, that dw/duy is small, we have the 
relationship 

oc n^^ . (2.9) 

That is, the neutrino mass increases as the neutrino density decreases. For w —1, the 
neutrino mass is simply inversely proportional to the neutrino density. 

2.1 Clustering and constraints 

In the previous discussion, we have made the two assumptions that a) the neutrino energy 
density is uniformly distributed, and b) we can find the value of the background field by 
minimizing its effective potential. Here we consider the validity of these assumptions, how 
deviations affect our scenario, and what additional constraints we need. 

We begin by considering the dynamics of the acceleron field A. So far we have assumed 
that we can treat the neutrino source term for A as being spatially constant. In order that A 
does not vary significantly on distances of the order of the inter-neutrino spacing, which, for 
thermal relic neutrinos, is currently 0(10^^ eV), the mass of the A particle at the present 
time must be less than ^ (10^'* eV). Such a light new particle could also couple to electrons 
and nucleons, and lead to an observable new force on scales longer than a millimeter — we 
consider expectations for such couplings and constraints in sec. ^. Another issue is whether 
the acceleron field remains at the minimum of its potential or oscillates. It is self-consistent 
to assume the acceleron field evolves adiabatically and remains at its minimum provided the 
time and distance scales over which the neutrino density varies are long compared with the 
oscillation period. A lower bound on the mass may be derived from the assumptions that the 
acceleron has evolved adiabatically since nucleosynthesis, and that the mass has either not 
changed significantly, or has increased since then, as is the case in the prototypical models 
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we discuss later. The acceleron mass should then be greater than the Hubble constant at 
nucleosynthesis, about 10"^"^ eV. There are various ways to explain how such a light field 
could have settled down to the minimum of its potential prior to nucleosynthesis [51] — as 
these are all quite model dependent we will not discuss them here^. 

We now consider the effects of neutrino inhomogeneity. Once neutrinos have become suf- 
ficiently nonrelativistic, they can cluster. Gravity will pull some neutrinos into existing dark 
matter halos. For instance, the authors of ref. [57] have studied the gravitational clustering of 
massive neutrinos in the background of dark matter halos modeled in existing high-resolution 
numerical studies, and found that significant overdensities can occur for massive neutrinos. 
Since in our framework, the neutrino mass is not a constant, one cannot directly appropriate 
the results. However, the vast majority of the clustering happens within a z of one, and a 
reasonable estimate is to use the neutrino mass at z = 1 as a benchmark. For instance, for 
a neutrino of mass 0.6 eV at z = 1, an overdensity of ~ 30 in the local group is reasonable. 
Since the cosmological neutrino density has changed by a factor of eight since that time, as- 
suming the neutrino mass is inversely proportional to the local density gives a neutrino mass 
in our vicinity today of ~ 0.6(8/30)eV ~ 0.15 eV. 

These numbers are rough estimates, but they illustrate an important point. The "cosmo- 
logical" neutrino mass (outside gravitationally bound systems) would then be 0{5 eV), which 
implies w « —0.8. At the same time, one would measure a mass on Earth of rrii, ~ 0.15 eV. 
Without considering the implications of clustering, such a low rrii, on Earth would seem in 
conflict with the equation of state in this model. Precise numerical simulations will be needed 
in order to determine the proper relationship between the local neutrino mass, the current 
cosmological neutrino mass, and cosmological data. 

A further possible source of neutrino clustering is the force between neutrinos mediated 
by the acceleron. These effects are quite model dependent, requiring knowledge of the range 
and strength of the force, and also numerical study, which is beyond the scope of the present 
paper. Here we simply note that clustering is suppressed by several factors, including the short 
range of the neutrino force and the nonnegligible neutrino velocity. In any case, for w near 
-1, the homogeneity of the dark energy is affected very little by neutrino clustering. Eqns. 

(l+w) 



and 2.£ show that the dark energy density varies with the neutrino density as nl, 



Neutrino clustering is of interest mainly due to its effects on neutrino mass phenomenology. 



^It is tempting to speculate that oscillations of the acceleron field could provide the dark matter. If 
the acceleron potential, averaged over an oscillation, is flatter than quadratic [53-55], acceleron dark matter 
would have negative pressure, similar to dark energy. However it would be unstable to clumping [53-56]. 
Such clumping would continue until the pressure vanishes, and subsequent evolution would be like normal 
pressureless matter. Also interesting to consider is the case of a steeper than quadratic potential, leading to 
positive pressure. One might hope [38,41,42,52] that obtaining dark matter and dark energy from the same 
source would help address the dark matter coincidence problem. However, the discussion of section 
that this is not the case. 



1.1 



shows 
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3. Models with a variable mass singlet 

A great deal of phenomenology can be learned directly from model independent features of the 
MaVaN scenario, but certain questions — specifically the dynamics — require a specific model. 
Dirac neutrinos require SM singlet fermions, and Majorana neutrino masses (a dimension five 
operator in the standard model), often arise from integrating out sterile fermions. A very 
simple implementation of this scenario is to allow the sterile neutrino mass to vary, causing 
the mass of the light eigenstate mass to change accordingly. 

As a general framework, we will assume that the active neutrino mass will arise from 
integrating out a heavier sterile state, whose mass depends on a light scalar field. We take a 
Lagrangian of the form 

C = mirl^ll^r + M (AjUrl^r + h.c. + log {I + \M (A)/ n\) , (3.1) 

where I'r is the sterile state, and will assume 

M{A)/n » 1 . (3.2) 

At scales below the right handed neutrino mass, we have the effective Lagrangian 

2 

^ = F(5) + + log(M(^)//x) . (3.3) 

In general, will find A ~ 10^'^ eV and naturalness constraints discussed later suggest m/^ ~ 
leV. 

From the discussion in section Q all models of this form will have the same equation of 
state, and the same neutrino mass for a give neutrino density. However, depending on the 
form of M{A) which we choose, the dynamics can be quite different. We will consider two 
limiting cases, one in which M{A) = XA and another where M(A) = Me-^^^-^^ . The scalar 
potential in the former case is very flat, while in the second it is just an ordinary quadratic 
potential. We will refer to these as the linear model and the exponential model, respectively. 
Notice that the second case is not unlike what might appear in some models where ^ is a 
dilaton field for a new force which is the dynamical origin of the mass term M{A). We leave 
consideration of the possible experimental signals of such a new force coupled to the sterile 
neutrino sector for future work. 

These two models should by no means be taken as the only possibilities. We choose 
them because they illustrate a wide range of possible dynamics which one might find in 
this framework. Although it is interesting to consider the the origin of these potentials and 
couplings, and possible connection with the other fermion masses, and other sectors of the 
theory, here we leave these issues aside. 

3.1 Quantum corrections and naturalness 

To have an equation of state near w = —1 one requires a rather flat potential for the neutrino 
mass. In general, quantum effects lead to corrections to the potential, which are proportional 
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Figure 2: One loop contribution to the rrii, potential. 



to the second and fourth powers of the masses of virtual particles. Thus, from quite model- 
independent considerations, we can see that the equation of state from theories of variable 
mass particles is strongly constrained by naturalness issues. For the moment, let us consider 
a variable mass particle with mass m, and the correction to the potential for m. 

At one-loop, we necessarily have the quadratically divergent contribution to the scalar 
potential in figure |2|. We will throw out the contribution proportional to A^, and assume that 
either the true cosmological constant is zero, or that IR dynamics of gravity eliminates its 
effects. We are left with a term 

^V„~|^, (3.4) 

where A is the cutoff of the theory. 

The dependence on A^ indicates that this correction is quadratically sensitive to the 
cutoff scale physics, and one might assume that this short distance physics for some reason 
is such that the contribution eqn. vanishes, although this is conventionally thought to be 
unnatural. However, there is also the cutoff insensitive correction to the potential 

'^^0 = -TjTTl-'^^logCW/^) : (3-5) 

where fi is the renormalization scale, which we should take to be of order the value of m 
today. 

If we want such a variable mass particle to give dark energy with an equation of state 
parameter w close to —1, we require 



V 



< 1 ^ < (lO-^eV)^ (3.6) 



As this argument applies to all particles of varying mass, we see heavy particles are 
generally unsuitable^ but that the neutrino is an ideal candidate for contributing to the 
dark energy. Conventional naturalness implies that the quadratic terms in the potential are 
sufficiently small only if the cutoff of the neutrino mass sector is less than of order an eV. 



^An exception would be if the heavy particle were a member of a nearly degenerate supermultiplet. 
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Examples of natural MaVaN models with a low cutoff for quadratic divergences include a 
neutrino mass due to integrating out a sterile neutrino which has mass near an eV and is 
either a member of some highly supersymmetric sector, or is an extended object, with a 
mass and compositeness scale of order an eV. In any event, naturalness suggests that new eV 
mass states would be detectable at the upcoming MiniBooNE experiment, which can test for 
neutrino oscillations into sterile states in this mass range. 

One might be concerned that this introduces another coincidence into the problem: 
namely that the present neutrino mass is fairly near the expected cutoff of the theory. How- 
ever, as total dark energy is set by the cutoff of the theory, this coincidence is equivalent 
to the fact that we live during the era when the neutrino mass is comparable to the fourth 
root of the dark energy, and equivalently to the accident that the cube root of the neutrino 
number density is within a factor of 10 of the fourth root of the dark energy. This is, in turn, 
directly related to the coincidence of the loosely similar values of the fourth root of the dark 
energy and the temperature of the microwave background, or to the fact that the radiation 
density was comparable to today's dark energy density, at a redshift of order 10. Since the 
neutrino mass is near the expected cutoff of the theory, we would expect that the present 
expansion will not last for long. Once the maximal value of the neutrino mass is reached, the 
scalar field is at its minimum and, since we have assumed vanishing cosmological constant, 
accelerated expansion will cease. The observational implications of this prophesy are limited 
to the expectation of new eV-scale states with couplings to the neutrino. 

4. Implications for Phenomenology, Cosmology and Astrophysics 

In section 2 we showed that if only nonrelativistic neutrinos of uniform number density are 
present, then their mass varies as . Here we consider the case where both background 
neutrinos and relativistic neutrinos are present. We will find situations where the extra 
neutrinos significantly change the mass from the background dominated case and examine 
the implications for neutrino mass limits. 

The contribution of neutrinos to the energy density is still given by eqn. ( |2.1| ) 



but now f{k) includes both background neutrinos and those from local sources. Since these 
sources are nuclear or high temperature "neutrinostrahlung" reactions, the additional neu- 
trinos are relativistic and we can write f{k) = fcuB{k) + frei{k)- The local neutrino mass is 
now given by 



The local neutrino mass will differ from its background dominated value when the second 
term is larger than or of the order of n'^'^^ . An estimate of this can be made by assuming 
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the relativistic neutrinos all have the same energy Ey = k ^ my. In this case the above 
expression becomes 

n^-^ + ^nr' + yo'(m.) = 0, (4.3) 

where n^y^ is the density of the relativistic neutrinos as seen in the background frame. 

Our approximate measure for finding a non-background mass is to take the ratio of the 
first two terms in (4.3), assuming that ruy does not change. For practical purposes we shall 
choose 1 eV as a value for rriy which may be an overestimate in some cases, but since we are 
seeking an approximate test for interesting environments, this is sufficient. Thus we define 

1 eV n:f 



(4.4) 



When this is large, the local neutrino mass differs from its value in a background dominated 
environment. Since Vq is a decreasing function of my, larger neutrino densities lower the 
local mass, and since the density cannot be lower than that of the background, neutrinos in 
a region dominated by the background have the heaviest mass possible. 

We will also take a number density of neutrinos Uy'^^ ~ 100 cm~^ = 8 x 10"^'^ eV^, which 
is the expected number for one species of thermally produced neutrino. As we have discussed, 
gravitational clustering can raise this value, which will, in turn, result in smaller values of r. 
Since we are merely interested in isolating those situations which are of possible interest, this 
is a conservative assumption. We shall use values of El'^^ and n^*^' which are characteristic of 
the environment under consideration. Note that the scaling properties of the neutrino mass 
already described are only applicable for non-relativistic neutrinos. If relativistic neutrinos 
dominate, the scaling will be considerably different. As this is particularly relevant in the 



early universe, we will discuss this in detail in section 4.2 



We will now use this measure to show that all neutrino environments on Earth are 
background dominated. Since r is an order of magnitude estimate, we will pass over the 
model dependent question of which neutrino fiavors couple to A, and use the flavor with the 
highest density to look for environments where relativistic neutrinos result in a lower mass. 

The most numerous solar neutrinos are the pp neutrinos, with energies of order 0.1 MeV 
and a flux on earth of 6 x lO^'' cm~-^s~^ [58]. Their number density is thus 2 cm^^ = 2 x 10^^"^ 
so r ~ 10 and the background neutrinos dominate solar neutrinos on earth. In the 
solar core r is larger, but we will discuss this in sections ( [4.1D and ([4.4| ). 

In reactor cores, the neutrino density can reach ~ 10^ cm~^ ~ 

10-11 eV^ [59], and the 
energies are again of order 0.1 MeV, so there again r is small (~ 10~^). 

In the most intense neutrino beams like those at K2K [60] and MiniBooNE [61], proton 
pulses collide with a target to produce pions that decay into GeV muon neutrinos. One can 
consider the limit in which a neutrino is produced for each proton on target. With roughly 
6 X 10^^ pot (protons on target) in 1.1 fis, with a source of cross section ~ 5 cm^, one yields 
an upper limit of 4 x 10'' v cm~^. For GeV neutrinos, this still yields an upper limit on r less 
than one. 
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Neutrinos from other sources such as the Earths radioactivity [62] and cosmic rays are too 
diffuse to be important. We note that nuclear fireballs might briefly produce neutrino densities 
large enough to change their mass, but no physics has been done with fireball neutrinos and 
we do not propose it should be. 

We conclude that the local background neutrinos dominate in all earth based experiments 
and that the tritium beta decay limit [63] gives the most stringent upper bound of ~ eV for 
the background neutrino mass. 

4.1 Interesting Environments 

For neutrino masses to vary significantly, we need either a high neutrino density relative to 
the background, or if few high energy neutrinos are present, then significant variation in the 



background itself. As discussed in section (2J), we will not consider background variations 



in depth, but we will now turn to several high density cases. They are in the early universe, 
supernovae, and stellar cores. 

Redshifting 100 v cm~^ back to the BBN era with a z of lO^'^, we find neutrino number 
densities of lO^'' eV^ and extremely light active neutrinos (r ~ 10^^). 

In core collapse supernovae, the early neutrino number density can reach 10'^^ cm~'^ (10^^ 
eV^) [64], which with order MeV energy gives an r of ~ 10^"^, certainly enough to lighten our 
active neutrinos. A separate issue from simply lowering the local neutrino mass is whether A 



emission changes the rate of cooling of the young neutron star, but we show in section (4.3) 
that the effect is small over broad regions of parameter space. 

While the above two environments are clearly dominated by relativistic neutrinos, stellar 
cores are a marginal case. The pp neutrinos that failed to dominate the background on 
earth are denser in the solar core. Within roughly 0.1 Rq [58], their density is ~ (1 AU/0.1 
Rq)"^ X 2 X 10~^^ eV^ = 9 X 10~^ eV^, and an energy of 0.1 MeV gives them an r of 
order 1. Higher energy neutrinos such as those from reactions have a negligible r due 
to both their lower number density and higher energy. Thermally produced neutrinos from 
pair reactions such as 'je — > euv have a lower energy but are also more diffuse and have a 
negligible r ~ 10~^ [65]. An r > 1 for any of the components will affect the masses of all 
neutrinos, of course, so the pp neutrinos may have an impact on the MSW properties in the 



sun. We will discuss this in greater detail in sections i.4 . 

Before discussing these cases further, we note two other environments in which neutrino 
masses could vary significantly. White dwarves cool largely by thermally produced neutrinos 
and can easily have neutrino densities of order 10 cm~ which gives a large r ~ 10^ (see 
Appendix C of [64]). Also, large stars such as red giants can produce more neutrinos from 
smaller regions than in the sun [64]. As the solar case is already marginal, we may find 
reduced neutrino masses in the cores of larger stars. However, since no observations have 
been made with neutrinos from these sources we will not treat them further here. 

Table 1 gives a summary of environments from which neutrino physics has been extracted, 
with r > 1 showing those in which the neutrino mass is lowered. 

We now consider the BBN, SN, and solar cases in slightly more detail. 
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Eu 


r 


BBN 


1030 


~ MeV 




SN cores 


1021 


~ MeV 


1027 


solar core 


~ 10-'^ 


0.1 MeV 


1 


reactor cores 


10-11 


0.1 MeV 


10-4 


beams 


< 3 X 10-7 


GeV 


< 10-2 



Table 1: r values for various environments, r > 1 implies that neutrino mass parameters could be 
modified. 



4.2 Limits on sterile neutrino and acceleron production from BBN 

Under ordinary circumstances, a sterile neutrino that has both a) appreciable mixing, and b) 
sufficiently small mass, will come into thermal equilibrium prior to nucleosynthesis. This can 
have serious consequences for BBN, making models with sterile neutrinos very constrained, 
cosmologically. Here, we have not only a sterile neutrino, but also the possibility of producing 
light A quanta, so one must take care that BBN is not radically changed. 

Concerning the sterile neutrino, the key point is that while parameters may be interesting 
today (i.e., accessible at the MiniBooNE experiment), these are effectively unrelated to the 
values at times prior to nucleosynthesis. It is important to note that the relationship rui, oc nj^ 
holds only for non-relativistic neutrinos. For relativistic neutrinos, we need to minimize the 
potential with the relativistic form for the scalar neutrino density, namely 

^ + Fo'(m.) = . (4.5) 

Following the same procedure as we did for non-relativistic neutrinos, we arrive at 

,2 



.,■11- "^^'^ - ^uVpi^u) . . 

3TV 3V ' ^ ' ' 

Now, assuming a slow variation in it;, we have 

oc a(3'"+i)/2 , (4.7) 

which is considerably slower than in the non-relativistic case. 

In the microscopic models presented, it was a sterile neutrino whose mass was varying, 
and hence, for a relativistic system, the sterile neutrino mass will vary differently in the hot 



early universe compared with later. In particular, if we consider the model in Eq. 3.1, then 



for relativistic neutrinos, the expectation value of for M{A) (the sterile neutrino mass) is 

MiA) ^ . (4.8) 

Since mir > rrii, > A the sterile state should be above the temperature scale and should not 
affect the light element abundances. 



- 13 - 



A 



A 
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Figure 3: Neutrino annihilation diagrams. 



We now turn our attention to the acceleron field which is setting the mass of the right 
handed neutrino. We certainly expect acceleron particles to be light (no heavier than 10~^ eV 
at the present era), and so they can be potentially thermalized in the early universe as well. 

Through the diagram of figure ^a, we expect an annihilation rate in the early universe 

4 

T , (4.9) 



1 



dA 



while through the diagram of figure we expect a rate 



1 

T 



(4.10) 



In addition, there can be A production through ^-strahlung following a weak interaction, 
which has a production rate 



'dA 



rp5 



(4.11) 



so the requirement that these particles are not in thermal equilibrium gives us 



drrii, 
~dA 



< 



< 



10 



-22 



10-34gV 



PI 



(4.12) 



'dA 



<1 



where the final relationship comes from evaluating these expressions at the era of BBN. 

Let us consider our two example models. In the linear model, we have during nucleosyn- 
thesis 



{A) 



AA2 



(4.13) 
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hence 



Hhr. = < 10-16 U 14) 



'dA 



< 10-22 eV 



hence neither contributions to the production of A will thermalize it. 
In the exponential model we have 



^ ^ ^ , (4.15) 



dA Tf 



dA^ T \ /4 

hence to avoid thermalization, one should have / ^ 450 ev. 

Interestingly, because the light neutrino mass eigenstate is becoming increasingly sterile 
at late times, the production of acceleron quanta can happen at late times. Should this occur 
when the neutrinos are non-relativistic, their number density would become exponentially 
suppressed, which would be problematic even for logarithmically flat potentials. Because the 
particles are Majorana fermions, there is an additional velocity suppression of their cross sec- 
tions when they are non-relativistic, and hence the annihilation rate of neutrinos drops faster 
than the Hubble rate. Hence, we need only check that they are not in thermal equilibrium at 
the scale when they go non-relativistic. In the models we are considering, this is when 

T^m„^ A^T (4.16) 

Since A ~ IQ-^eV, this occurred at z ~ 10 at which time the Hubble scale was different by a 
factor of roughly ten. 

We begin by considering the linear model. At the non-relativistic transition, the annihi- 
lation rate goes as from both diagrams goes as 

1 A^A^ 



(4.17) 



And requiring this to be less than the Hubble scale requires 



^>4eV (4.18) 

So for A ~ 1/100 the scenarios discussed should be safe. For larger values of A and lower 
values of mij- a more rigorous analysis should be performed. 

In the exponential model, at the transition, we have (taking [A) ^ f, which is natural if 
the sterile neutrino mass is to be of order eV.) 

1 A^ 

- ~ - (4.19) 
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Requiring this to be less than Hubble at 2; ~ 10, we yield 



/ > 10 keV , 



(4.20) 



which is more restrictive than previous limits. 
4.3 Supernovae 

Supernovae are interesting to us for several reasons. In contrast to neutrinos from white 
dwarves and red giants, neutrinos from a core collapse supernova have been observed [66-68], 
and we must check that A emission does not cool a protoneutron star too quickly to alter 
the observed spectrum of neutrinos emitted in the first 10 seconds. Further, if r-process 
nucleosynthesis occurs in the envelopes of core collapse supernovae, it will be sensitive to the 
flavor composition of the neutrino wind emanating from the newly formed neutron star. 

We can quickly check that A emission is not significant. Neutrinos inside the early 
neutrino sphere have energies of order 10 MeV, and we need emission rates of less than 1/10 
sec~^. This requires only slightly stronger limits as those from BBN: 



These are still easily satisfied by the linear model. The v4-strahlung limit is itself tighter 
than the BBN limit, but it is again redundant given the limit from uv AA. Thus A 
emission does not speed the early cooling of supernovae. 

A cooler neutrino sphere persists for roughly 10-100 years [64], during which time neutrino 
cooling dominates in the standard picture. This again may be altered by A emission in our 
case, and the longer timcscalc sets somewhat tighter limits than those imposed in the first 
10 seconds. However, since cooling rates of young neutron stars have not been accurately 
measured, due to the difficulty in separating X-rays from the stars surface with other sources 
[69], we will not pursue this further here. 

The r-process is thought to occur in the envelopes surrounding core collapse supernovae, 
where neutron capture onto heavy nuclei proceeds more rapidly than the beta decays of these 
nuclei. The process is sensitive to the flux from the core, since for example the fg's can turn 
the free neutrons into protons, reducing the neutrons available for the r-process. Since these 
protons end up inside a particles this problem is known as the "alpha effect" [70,71]. Both 
active-active oscillations [72, 73] and active-sterile oscillations [74-76] have been suggested as 
important mechanisms for both the r-process and shock dynamics, and it is possible that 
active-sterile mixing with a dm^ of 0.1 - 100 eV^ and matter enhancement might solve the 
above "alpha problem" [77]. 




(4.21) 
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In our variable mass singlet model, we have a sterile state whose mass would be ~ 
MeV in the core, that becomes progressively lighter as it propagates outwards into the less 
dense envelope, and an active state that get progressively heavier and more sterile. Should 
MiniBooNE see a positive result, its implications within the MaVaN scenario for the r-process 
could be significantly different and would need to be carefully analyzed. 

4.4 Solar neutrinos 

Our estimate for pp neutrinos in the solar core shows that a detailed calculation is needed in 
the event that A couples to Ue- If the electron neutrino mass is modified, this could affect our 
interpretation of the KamLAND and SNO results^. For neutrinos with constant mass, fgS 
are produced in the sun's core and oscillate into Ufj^s and i^t-s via the MSW effect whilst still 
inside the sun. In our scenario, background neutrinos dominate the acceleron potential until 
we are sufficiently near the core of the sun, at which point the mass of the electron neutrino 
could be lowered depending on it's energy. There is the possibility that the neutrino mass 
matrix would be altered near the solar core, modifying the interpretation of the solar neutrino 
results, but this is model dependent and the detailed consequences are yet to be worked out. 

4.5 Baryogenesis 

There are many baryogenesis scenarios which constrain neutrino masses — see for example refs. 
[78-85]. Majorana neutrino masses violate lepton number (L), and, since high temperature 
nonperturbative electroweak processes violate baryon number (B) but conserve B-L, an upper 
bound on Majorana neutrino masses, which depends on the baryogenesis scale, can be derived. 
An attractive baryogenesis mechanism which constrains both the sterile and the light neutrino 
masses is to use CP violation in the decays of heavy sterile neutrinos to create L, which 
is subsequently partially converted in B by nonperturbative electroweak interactions. The 
constraints obtained on neutrino masses from baryogenesis scenarios all assume constant 
neutrino mass and will clearly have to be reevaluated in the MaVaN scenario. In the models 
of sec. ^, early on the sterile neutrinos were much heavier and the active neutrinos much 
lighter, and so the constraints from baryon number washout on neutrino masses become very 
weak. 

5. Other Roles for the Acceleron? 
5.1 Long Range Forces 

Could the acceleron mediate forces between other particles than neutrinos? There are strong 
limits on new forces for between nucleons and/or electrons at scales longer than 100 microns 
[51]. For instance, the operators 

me{A)ee , md{A)dd , mu{A)uu , -oTT^G" G'''^"' . (5.1) 

similar effect has been studied in the context of "neutrino clouds" [47], although the effect is due to 
clustered primordial neutrinos, not pp neutrinos. 
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could lead to violations of the force law and of the weak equivalence principle, on scales 
longer than a millimeter, if m'^,m'^,m'^, or g'^ ^ 0. The Yukawa force mediated by, e.g., the 
operator me{A)ee , is proportional to m-g^, and has range 

Standard Model quantum corrections will necessarily induce acceleron dependence into 
such operators, but with negligible coefficients, suppressed by loop factors and by a factor at 
least as small as 

. 

In the explicit models of sec. ^ the minimum natural size of any such acceleron mediated force 
is too small to give observable deviations from the inverse square law or the weak equivalence 
principle^. 

Although we lack understanding of the role of A in the fundamental theory, we can 
postulate that it could be a modulus for other parameters of the standard model besides 
neutrino masses, and test for this by searching for violations of the weak equivalence principle 
and of the inverse square law at distance scales longer than a millimeter. The observation or 
nonobservation of such effects constrain the mass and couplings of the acceleron, and guide 
model building towards its fundamental origin. 

5.2 Varying a 

Besides leading to new forces, since the value of A depends on s,y, nontrivial A dependence of 
the operators implies that "constants" of the standard model depend nontrivially on the 
scalar neutrino density. This is different from typical treatments of varying constants, which 
mainly consider temporal variations (for reviews see ref . [88-90] ) . It is particularly interesting 
to consider the operator 

' (5-3) 

in light of the recent astrophysical evidence from quasar absorption lines that a may be 
variable [91,92]. 

Such dependence is highly constrained by BBN, since Sjy was larger then by a factor of 
at least 10^^ . The success of BBN implies that the constants of nature have changed little 
since. One might think that the BBN implies that no interesting observable dependence of a 
or other constants on Siy is possible. However, if the dependence of a on Sj/ is sufficiently flat, 
e.g. logarithmic, then, given that Si^ varies steeply with redshift and is likely to vary with, 
e.g. the matter density due to gravitational neutrino clustering, it is possible to reconcile the 
observation of astrophysically varying a with BBN. 

The hypothesis that a depends on explains a troubling aspect of the recent claims. 
It is calculated that the presence 2 Gyr ago of a natural nuclear reactor on Earth implies 



^Although not interesting for long range forces, the coupling of the acceleron to the electron which is 
induced by quantum corrections can be large enough to induce phenomenologically important changes in the 
masses of neutrinos propagating through electron dense matter such as the core of the sun or solids on the 
earth. This effect differs from the usual matter effect [86] in being the same for neutrinos and anti neutrinos 
and in not being neutrino energy dependent. A phenomenological study of such effect is underway [87]. 
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that a cannot have fractionally changed by more than ~ 10~^ since then [93]. Ref. [91] 
interprets its data as evidence that a has increased by a fraction ~ 10~^ since redshift of 
greater than ~ 0.5. If one assumes linear variation of alpha with time, the results of refs. [91] 
and [93] are inconsistent. However, the local value of should have changed little since the 
formation of the galaxy, but could be quite different outside our galaxy or galaxy cluster, 
and so the terrestrial value of a could be essentially constant while the value of a outside the 
galaxy could be slightly varying. To test this, one could look for dependence of a on, e.g.^ 
the matter density, rather than on redshift. Similar effects have been studied in some mixed 
quintessence/ varying a models [94,95]. 

5.3 MaVaNs and High Energy Cosmic Rays from the Z-Burst Scenario 

One of the most puzzling questions in high-energy astrophysics is the nature of the highest 
energy cosmic rays [96-103]. The expected Greisen-Zatsepin-Kuzmin (GZK) cutoff [104,105] 
in their energy has not been observed [106], suggesting the possibility of an exotic origin of 
the highest energy cosmic rays. 

One possibility is the Z-Burst scenario [107] in which ultra-high energy neutrinos reso- 
nantly annihilate with thermal relic neutrinos. The possible success of this model depends 
strongly on the neutrino mass, both to yield cosmic rays of the proper energy, and to in- 
duce sufficient clustering to amplify the sources. However, the 68 % preferred mass range of 
0.08 eV- 1.30 eV (extragalactic source for UHE z/'s) or 2.1 eV-6.7 eV (halo source) [108] is 
significantly different from the measured atmospheric neutrino splitting. If all the neutrinos 
are degenerate, much of the parameter space is in conflict with limits on the neutrino mass 
from structure formation. 

In the MaVaN scenario, there are two significant differences from ordinary neutrinos^. 
The first is what has already been stated: the cosmological limits on neutrino masses simply 
do not apply here. Hence the entire parameter space is opened up. The second, arises from 
the different density of neutrinos in gravitationally bound systems. As we have stated the 
measured value of neutrino mass on the Earth need not be the same as in other systems, nor 
does the mass in the interior of a cluster need to equal the mass near the exterior. This would 
in principle require an integration over a continuum of mass values, rather than a particular 
value, as in the fixed mass theory. As a consequence, the parameter space for the Z-burst 
scenario as it relates to terrestrial measurements of neutrino mass is opened considerably. 
To properly understand this, a careful analysis of neutrino clustering along the lines of [57] 
would be essential. 

6. Outlook and Smoking Guns 

Inspired by the similarity between the the dark energy and neutrino mass scales, we have 

®The effect of neutrino clustering has also been studied in the context of neutrino clouds [48] where consid- 
erable changes to the Z-burst scenario also occur [109]. The physics is quite different, however, as the changes 
are essentially due to the large Fermi momenta of the neutrinos inside the neutrino cloud. 



- 19 - 



proposed a link between the dark energy and the neutrino mass, which imphes that the 
neutrino mass is approximately inversely proportional to the scalar neutrino density. The 
neutrino mass varies in such a way that the total energy density in cosmic neutrinos is always 
similar to the dark energy density. The fact that the neutrino mass today is not very different 
from the fourth root on the dark energy is related to the coincidence that we live during the 
epoch when the microwave background temperature is within an order of magnitude of the 
fourth root of the dark energy. 

To test the scenario terrestrially requires an extremely intense source of neutrinos — 
intense enough so that the scalar density, which is suppressed by niu/Eu, is larger than the 
scalar density of the cosmological background neutrinos. To date, the only sufficiently intense 
terrestrial source is a nuclear bomb, in which we see no way to test for neutrino masses. There 
are however several possible future observations which could provide distinctive evidence for 
MaVaNs. One possibility is that the upcoming MiniBooNE experiment [22] could provide ev- 
idence for one or more sterile neutrinos. Since, assuming fixed neutrino masses, the mass and 
mixings of such detectable sterile neutrinos would likely be problematic for nucleosynthesis, 
this alone would suggest that the neutrino mass should vary. For solutions to the questions 
relating to sterile neutrinos in BBN, in particular those involving primordial lepton number or 
additional fields see [84,110-116]. Another possibility is that cosmological or supernova con- 
straints on neutrino masses could eventually conflict with the neutrino masses derived from 
observations within our solar system-again implying that the neutrino masses should vary. 
Finally, an acceleron coupling to the electron or nucleon could lead to an observable matter 
density dependence of neutrino oscillations which differs from the usual matter effect [87]. 

There are two model independent predictions of MaVaN dark energy. One prediction 



is the relation of eq. 2.7 between the dark energy equation of state parameter and the cos- 
mological neutrino mass. To test this, one needs to precisely measure the equation of state 
parameter. This will determine the cosmological neutrino mass today. To obtain a prediction 
for the local neutrino mass, which depends on the neutrino density in our vicinity, one needs 
to use detailed numerical modelling to study neutrino clustering in our galaxy. The second 
prediction is that if all neutrino masses vary in the same way, neutrino masses at high red- 
shift scale like (1 + z)^"', and so high redshift cosmological data should show no evidence for 
neutrino mass. 

Still requiring much further study are the consequences of MaVaNs for neutrino oscillation 
experiments, solar neutrino physics, supernovae, neutrino clustering, structure formation, 
baryogenesis, and high energy cosmic rays. Also interesting is the cosmology of and the 
origin of the acceleron, the scalar field which mediates the dark force between neutrinos. 
Searches for deviations from the equivalence principle and the 1/r^ force law at scales of 
order a millimeter or longer could shed some light on the acceleron, as could variation of 
other 'constants' with the neutrino density. 
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